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SULPHUR CHEMISTRY IN THE ENVELOPES OF MASSIVE YOUNG STARS 
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Abstract 

We present submillimeter observations of SO, SO2, 
H 2 S, H 2 CS, OCS, NS and HCS+ toward nine massive 
young stars. The outflow contributes «50% to the SO 
and SO2 emission in 15 — 20" beams, more than for CS, 
where it is 10%. The SO2 abundance increases from dark 
cloud levels in the outer envelope (T < 100 K) to lev- 
els seen in hot cores and shocks in the inner envelope 
(T > 100 K). Molecular abundances are consistent with 
a model of ice evaporation in an envelope with gradi- 
ents in temperature and density for a chemical age of 
~ 30000 yr. The high observed abundance of OCS, the 
fact that r ox (OCS)>T ex (H 2 S), and the data on solid 
OCS and H 2 S all suggest that the major sulphur carrier 
in grain mantles is OCS rather than H 2 S. For most other 
sulphur-bearing molecules, the source-to-source abundance 
variations by factors of up to 10 do not correlate with 
previously established evolutionary trends in temperature 
tracers. These species probe the chemically inactive outer 
envelope. Our data set does not constrain the abundances 
of H 2 S and SO in the inner envelope, which, together with 
S0 2 , are required to use sulphur as a clock. 
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1. Introduction 

Spectral line surveys at submillimeter wavelengths have 
revealed considerable chemical differences between star- 



forming regions (see van Dishoeck 2001 for an overview). 
While these differences indicate activity, the dependence 
of molecular abundances on evolutionary state and phys- 
ical parameters is poorly understood. Better insight into 
this relation would be valuable for probing the earliest, 
deeply embedded phases of star formation where diagnos- 
tics at optical and near-infrared wavelengths are unavail- 
able. This is especially true for the formation of high-mass 
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Figure 1. Examples of the spectra taken with the JC'MT. 
The top and bottom frequency scales of the left panel are 
the two receiver sidebands; the right panel only shows 
lines from one sideband. Adopted from van der Tak et 
al. (2002). 



stars, for which the order in which phenomena occur is 
much less well understood than in the low-mass case, and 
for which the embedded phase is a significant fraction 
(»10%) of the total lifetime of -10 6 yr. 

Sulphur-bearing molecules are attractive as candidate 
tracers of early protostellar evolution (Hatchell et al 1998). 
Models of "hot core" ch emistry (pharnlcy 1997 ), as well 
as of shock chemistry ( Pineau des Forcts et al 1993j ) 
predict strong variations in the abundances of sulphur 
molecules on time scales of ~10 4 yr. To test these sce- 
narios, this paper discusses submillimeter observations 
of sulphur-bearing molecules towards nine regions at the 
earliest stages of high- mass star formation, with luminosi- 
ties 1 x 10 4 - 2 x 10 5 L© at distances of 1 - 4 kpc. This 
paper only summarizes our results; a complete descrip- 



tion is given by van der Tak et al (2002) 



Proceedings of the conference "Chemistry as a Diagnostic of Star Formation, " University of Waterloo, Canada, 21-23 
August 2002 (C.L. Curry & M. Fich eds.). 
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Figi re 2. Velocity profiles of the SO 6e -5s line in the 
sources AFGL 2591 and W3 IRS5. 
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Distributions of the SO and SO2 abundances 
2591 for chemical ages of 300, 3000, 30000 
yr (curves from bottom to top) after Doty 
. The t ~ 3 x 10 4 yr model fits our data best. 



Spectra in the 230 and 345 GHz bands were taken in 
1995-1998 with the JCMT. The overview spectra in Fig. [I] 
show source-to-source differences in the relative and in 
the absolute strengths of lines of sulphur molecules. The 
line profiles (Fig. |^) reveal the presence of a narrow com- 
ponent, due to the envelope, and a broad one, due to the 
outflow. For SO and SO2, the outflow contributes «50% 
of the line flux. This fraction is much larger than for CS, 
where it is only 10%. Lines of other molecules are too 
weak to separate the two components. 

3. Results 

To calculate molecular abundances, we have used the 



(2000), and models of the temperature and density struc- 
ture of the sources (van der Tak et al 2000[ ). Table [l] 
compares our results with abundances in other types of 
sources. The best overall match is with "hot cores" . 

For some sources, the detection of many SO2 lines 
from a wide range of energy levels enables us to derive 
radial abundance profiles. Models with a constant SO2 
abundance of 10~ 9 fit the lines arising < 100 K above 
ground. However, the higher-excitation lines are under- 
produced, which indicates a 'jump' in the SO2 abundance 
by a factor of ~100 at a temperature of ~ 100 K. Infrared 



Monte Carlo program by Hogerheijde & van der Tak 



observations of SO2 absorption by Keane et al (2001) also 
indicate abundances of ^10~ 7 for the inner envelopes. 
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Table 1. Molecular abundances. 



Source 


CS H 2 


CS 


H 2 S 


HCS+ 


NS 


OCS 


SO 


S0 2 




10 -9 10 


-10 


10" 9 


10 -io 


10" 11 


10" 9 


lO" 9 


10" 9 


t i to TT> o tr 

W3 1RS5 


5 


3 


9 


0.2 




0.5 


10 


10 


W33A 


5 


13 


8 


2 


2 


20 


1 


1 


AFGL 2136 


4 


3 


<2 


1 


<0.5 




1 


0.5 


AFGL 2591 


10 


3 


20 


2 


0.5 


10 


10 


2 


S140 IRS1 


5 


7 


4 


5 


1 


<2 


2 


1 


NGC 7538 IRS1 


10 


10 


8 


2 


1 


2 


5 


1 


NGC 7538 IRS9 


10 


10 


8 


5 


1 


2 


1 


0.5 


MonR2 IRS3 


5 


7 


<4 








5 


1 


NGC 6334 IRS1 


10 


7 


10 


2 


10 


50 


2 


2 


Hot cores" 


8 


6 


9 




80 


5 


4 


20 


PDRs'' 


20 




6 








9 


0.1 


Dark clouds c 


1 


6 


0.8 


0.6 




2 


20 


4 


Shocks'* 


4 


8 


4000 


0.1 




10 


200 


100 


a Source average from Hatchell et al (1998) b Orion Bar: Jansen et al (1995) 


L 134N: Ohishi et al (1992) T 


Orion Plateau: Button et al (1995); Minh et al (1990) 



4. Chemical implications 

Most hot core models only consider single values for the 
temperature and density, while for our sources, these pa- 
rameters vary stro ngly along the lin e of sight. Using the 
UMIST database, |Doty et al (2002)| calculated the time- 
dependent gas-phase chemistry for the temperature and 
density structure of AFGL 2591. The model has most 
sulphur initially in S at T < 100 K and in H 2 S at T > 
100 K, to mimic the effect of ice evaporation. The re- 
sulting 'jump' in the abundances at r = 2 x 10 16 cm is 
consistent with our observations of S0 2 . See Fig. 3 and 
Doty (this volume). 

The model reproduces most of our observed abun- 
dances, and also the increase of the SO2 abundance at 
T> 100 K. However, OCS is underproduced by factors of 
10-100, suggesting that most of the sulphur in the grain 
mantles is in the form of OCS. This idea is supported by 
the high excitation temperature of OCS (100 K), much 
higher than the 25 K found for H2S. In addition, solid 
OCS has been detected in the infrared towards W33A, 
while H2S has not. 

The current data set is not suitable to trace chemical 
evolution within our sample. The CS, H 2 CS, HCS + and 
NS molecules mainly probe the chemically inactive outer 
envelope (T < 100 K). While OCS, SO and S0 2 show 
larger source-to-source variations, only for S0 2 , the data 
allow to determine the abundance in the inner envelope 
(T < 100 K) where ices are evaporating. The situation for 



H 2 S is unclear: models predict strong abundance varia- 
tions, but the present data do not trace the warm gas. 
Future observations of high-excitation lines of H 2 S, SO 
and OCS could exploit the full potential of sulphur as 
chemical clock. 
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